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I Abstract. This talk is a brief report on our results in high parton density QCD. The 

. main goal of this talk is to discuss the new evolution equation which was obtained 

I and solved in Refs. [1] [2] and to share with you the physics can be recovered by this 

00 ' equation. 

INTRODUCTION 

Oh! 

^ . In this talk, we briefly outline the main ideas, approaches and results which 

we developed in high parton density QCD (hdQCD). High parton density QCD 
involves a system of partons at short distances in which the QCD coupling constant 
as is small, but the density of partons becomes so large that we cannot apply the 
usual methods of perturbative QCD ( pQCD ). In essence, the theoretical problems 
here are non perturbative, but the origin of the non perturbative effects does not 
lie in long distances and large as-, which are present in the confinement region. 
First, we need to answer the question: 
1. Which parton density is high? 

The quantative estimates of which density is high, we can obtain from the s- 
channel unitarity [3] which we can use in two different form: 

1- o'tot{'y*p) < «em2 7ri?^, where R is the size of the target, aem is the fine 
structure constant ( see Ref . [6] for explanation why we have aem ) • This constraint 
can be rewritten as [3] [4] [5] 

3TcasxG{x,Q^) 
" = Q^^ ^ ' ■ 




2. cr^^iYp) ^ (^tot{l*p), where DD stands for diffractive dissociation. This 
inequaUty gives Eq. (1). 
Therefore, 

1. if K is very small <^ 1 ), we have a low density QCD in which the parton 
cascade can be perfectly described by the DGLAP evolution equations [7] ; 

2. if K < 1, we are in the transition region between low and high density QCD . 
In this region we can still use pQCD, but have to take into account the interaction 
between partons inside the partons cascade; 

3. if K > 1, we reach the region of high parton density QCD. which we are 
going to discuss here. 

Before doing this we want to make a remark on parton densities in a nu- 
cleus. Taking into account that for nucleus Ra = Rn x and xG Aix^Q"^) — 
AxGNiXjQ"^), we can rewrite Eq. (1) as 

37rasAxGN{x,Q^) a\ ^ ^ ro\ 
ka — 2 — A3 kn < 1 ■ (2) 

Therefore, for the case of an interaction with nucleii, we can reach a hdQCD region 
at smaller parton density than in a nucleon (see Fig.l). 

Fig.l gives the kinematic plot {x, Q"^) with the line which shows that the 

hdQCD effect should be seen at HERA. 

2. Two different theoretical approaches 

To understand these two approaches we have to look at the picture of a high 
energy interaction in the parton model. In the parton approach the fast hadron 
decays into point-like particles ( partons ) long before ( typical time r oc ^ ) 



FIGURE 2. The high energy interaction in the parton model. 



the interaction with the target. However, during this time r, all partons are in the 
coherent state which can be described by means of a wave function. The interaction 
of the slowest ( "wee" ) parton with the target completely destroys the coherence of 
the partonic wave function. 

The total cross section of such an interaction is equal to 

CTtot ^ N X ao (3) 

where 

m N = flux ( renormalized ?! ) of "wee" partons ; 

• (Jo = the cross section of the interaction of one "wee" parton with the target . 

One can see directly from Fig. 2 that the number of "wee" partons is rather large 
and it is equal to 

N oc e^"^ = with < n > ^ cuq \n{l/x) ; cuq = Cas ■ (4) 

We have to renormalize the flux of "wee" partons, since the total cross section is 
the number of interactions and if one has several "wee" partons with the same 
momenta they only give rise to one interaction. 

UN 1, we expect that the renormalization of the flux will be small, and we 
use an approach with the following typical ingredients: 

• Parton Approach; • Shadowing Corrections; • Glauber Approach; • Reggeon- 
like Technique; • AGK cutting rules . 

However, when N ^ 1, we have to change our approach completely from 
the parton cascade to one based on semiclassical field approach, since due to the 
uncertainty principle ANA(f) ^ 1, we can consider the phase as a small parameter. 
Therefore, in this kinematic region our magic words are: 

• Semi-classical gluon fields; • Wiezsacker-WiUiams approximation; • Effective 
Lagrangian for hdQCD; • Renormalization Wilson group Approach. 

It is clear, that for N ^ 1 the most natural way is to approach the hdQCD 
looking for corrections to the perturbative parton cascade. In this approach the 
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FIGURE 3. The current situation in the battle for hdQCD. 



pQCD evolution has been naturally included, and it aims to describe the transition 
region. The key problem is to penetrate into the hdQCD region where k is large. 
Let us call this approach "pQCD motivated approach " . 

For N ^ 1, the most natural way of doing is to use the effective Lagrangian 
approach, and remarkable progress has been achieved both in writing of the explicit 
form of this effective Lagrangian, and in understanding physics behind it [8]. The 
key problem for this approach was to find a correspondence with pQCD. This 
problem has been solved [9]. 

The hottest news of this conference is that these two approaches give the same 
evolution equation for the gluon structure function as reported by A. Kovner [10]. 

Fig. 3 shows the current situation on the frontier line in the offensive on hdQCD. 

3. Pluses and Minuses of pQCD motivated approach 

The goal of this talk is to show you how we derive this evolution equation and 
explain what physics we can recover using this approach. 

Let me list here the main advantages and shortcomings of this approach. I believe 
that the honest report on them will lead us to a deeper understanding of the main 
theoretical problems in hdQCD. 



1. Clear physical meaning of everything what we are doing; 

2. Transparent parameters that we are using for a selection of the Feyman 
diagrams ; 

3. Obvious interrelation with usual DGLAP evolution and the Reggeon Calculus 
( "soft" phenomenology ) ; 




4. Understandable limitation of the approach which actually means that we 
guarantee the accuracy of our theoretical calculations. In other word, this approach 
has reached a stage of a normal theory with a regular procedure of how to calculate 
corrections with respect to small parameters. 

5. It has a clear feedback for such interesting problems as higher twist contribu- 
tions, violation of factorization and so on. 

There is no hope to obtain any formal ( operator ) approach to the problem and, 
therefore, to find a formal justification of the resummation technique. It means, 
that we can prove everything only in the framework of the Feyman diagrams and in 
the framework of the operation rules with them. However, it means that we ignore, 
at least partly, non perturbative effects and that we are not able to formulate a 
procedure of calculation of such effects. 

Let me illustrate why it is very important an operator proof giving two examples. 

1. In the pQCD motivated approach we heavily use the so called leading log 
(1/x) approximation in pQCD ( LL(x)A ). The operation rule for LL(x)A is: 

xG'"W^(x,g2) ^ ZimU^oE^=i a§M„(x,g2) (5) 
= E^=i a^lim\,-.oMr,{x,Q^) = E^=i C„(Q2) ( ln(l/x) ) V 

We used the gluon structure function as an example. Eq. (5) shows that we changed 
the order of two operation : summation over n and limit at a; — > 0. Do we have a 
proof for this? As far as I know, the answer is no. We consider as a strong argument 
in favor, the fact that the calculation in the next-to-leading order in LL(x)A gives 
a reasonable answer. The real proof has been done only recently [9] in the effective 
Lagrangian approach, which leads to the LL(x)A of pQCD. 

2. One of the main non perturbative observable which enters any pQCD moti- 
vated calculation is so called the size of the target R ( see Eq. (1) ). We know the 
meaning of this parameter, namely, this is a correlation radius between two gluon 
inside a target at x ~ 1. In the effective Lagrangian approach we will be able 
to find an operator expression for this number, and to ask our friends from lattice 
QCD to calculate this value. I think this is the only reasonable way of how to take 
non perturbative QCD into account, at least, at present time. 



NEW EQUATION 

4- The GLR equation 

The first attempt to take into account the parton - parton interaction in the 
pQCD motivated approach was done long ago [3]. It was based on the simple idea 
that there are two processes in a parton cascade (see Fig. 2) : (i) the probability 
of the emission of an extra gluon is proportional to ag p where p is the density of 
gluon in the transverse plane, namely 



and (ii) the annihilation of a gluon, or in other words a process in which the 
probabihty is proportional to x a annihilation- cr annihilation can be estimated as 

(J annihilation oc o^r^, wheie T is the size of the parton produced in the annihilation 
process. For deep inelastic scattering, oc Therefore, in the parton cascade 
we have 

Emission (1 — > 2 ) : probability = P'^™''*'"'^^ oc as P \ (7) 
Annihilation (2^1): probability = panmhUation ^ ^2 ^ ^2 _ ^g-) 

At a; ~ 1 only emission of new partons is essential, because p <^ 1 and this 

emission is described by the DGLAP evolution equations. However, at a; ^ the 

value of p becomes so large that the annihilation of partons becomes important, 

and so the value of p is diminished. The competition of these two processes we 

can write as an equation for the number of partons in a phase space celP ( Ay = 
A ln(l/x) A ln(gVg2) ). 

^ P 7^ P ' (9) 



'd\n{l/x)d\n{QyQl) tt ^ Q 

or in terms of the gluon structure function xG{x, Q"^) 

d^xGix, Q^) as N. ^^^^^ _ ^ ^ ^^^^^ _ ^^^^ 



91n(l/,T)ain(gVg2) ^ v >^ ttR^ q^^ 

This is the GLR equation which gave the first theoretical basis for the consideration 
of hdQCD. This equation describes the transition region at very large values of Q^, 
but a glance at Fig.l shows that we need a tool to penetrate the kinematic region 
of moderate and even small Q^. 
5. Glauber - Mueller Approach 

We found that it is very instructive to start with the Glauber approach to SC. The 
idea of how to write the Glauber formula in QCD was originally formulated in two 
papers Ref. [11] and Ref. [5]. However, the key paper for our problem is the second 
paper of A. Mueller, who considered the Glauber approach for the gluon structure 
function. The key observation is that the fraction of energy, and the transverse 
coordinates of the fast partons can be considered as frozen, during the high energy 
interaction with the target [11] [5]. Therefore, the cross section of the absorption 
of gluon(G*) with virtuality and Bjorken x can be written in the form: 

<tiG*) = (11) 

I'dz l^n'{Q',n,x,z)a^{x,r^) [*f (g^ r„ x, .)]% 

where z is the fraction of energy which is carried by the gluon, is the wave 
function of the transverse polarized gluon, and cxNix, rf) is the cross section of the 



^ We will argue why we chose this cell below in section 6. 



interaction of the GG- pair with transverse separation with the nucleus. Mueller 
showed that Eq. (11) can be reduced to 



Jx X J -Xy TlTf JO TT ^ ^ 



with 



= r^G^«-^^(x, i) 



_ 6? 

and profile function S{bt) chosen in the Gaussian form: S{b^) = e «^ . 

Obviously, the Mueller formula has a defect, namely, only the fastest partons 
(GG pairs) interact with the target. This assumption is an artifact of the Glauber 
approach, which looks strange in the parton picture of the interaction. Indeed, in 
the parton model we expect that all partons not only the fastest ones should interact 
with the target. At first sight we can solve this problem by iteration Eq. (12) ( see 
Ref. [1] ). It means that the first iteration will take into account that not only the 
fastest parton, but the next one will interact with the target, and so on. 

6. Why equation? 

We would like to suggest a new approach based on the evolution equation to 
sum all SC. However, we want to argue first, why an equation is better than any 
iteration procedure. To illustrate this point of view, let us differentiate the Mueller 
formula with respect to y = ln(l/a;) and ^ = InQ^. It is easy to see that this 
derivative is equal to 

^ ±ldbni- e-<^(--^)^(^?) } . (13) 

The advantages of Eq. (13) are 

• Everything enters at small distances ; 

• Everything is under theoretical control ; 

• Everything that is not known ( mostly non perturbative ) is hidden in the 
initial and boundary condition. 

Of course, we cannot get rid of our problems changing the procedure of solution. 
Indeed, the non-perturbative effects coming from the large distances are still impor- 
tant, but they are absorbed in the boundary and initial conditions of the equation. 
Therefore, an equation is a good ( correct ) way to separate what we know ( small 
distance contribution) from what we don't ( large distance contribution). 

7. Equation 

We suggest the following way to take into account the interaction of all partons 

in a parton cascade with the target. Let us differentiate the fe^-integrated Mueller 
formula of Eq. (12) in y = ln(l/a;) and ^ = ln((5^/(5o)- It gives 



9^^!^ = { C + Hkg{x',Q')) + E,{kg{x',Q% , (14) 



where k.g{x, Q"^) is given by 



<^,Q') = ^^2^G{x,Q'). (15) 

Therefore, we consider on the l.h.s. of Eq. (14) as the observable which is 
written through the solution of Eq. (14). 

Eq. (14) can be rewritten in the form ( for fixed as ) 

— ^-^7 \ a ^ + ln(/«G) + Ei{kg)} = F{kg) ■ (16) 

oyo^ ay tt 

This is the equation which we propose. 

8. Nice properties of the equation: 

This equation has the following desirable properties: 

1. It sums all contributions of the order {asUiY) absorbing them in xG{y,^), 
as well as all contributions of the order of /t". Therefore, this equation solves 
the old problem, formulated in Ref. [3], and for Nc ^ oo Eq. (16) gives the 
complete solution to our problem, summing all SC ; 

2. The solution of this equation matches with the solution of the DGLAP evolu- 
tion equation in the DLA of perturbative QCD at k ^ ; 

3. At small values of k ( k < 1 ) Eq. (16) gives the GLR equation. Indeed, 
for small k we can expand the r.h.s of Eq. (16) keeping only the second term. 

Rewriting the equation through the gluon structure function wc obtain the 
GLR equation [3] with the coefficient in front of the second term as calculated 
by Mueller and Qiu [4] ; 

4. The first iteration of this equation gives the Mueller formula ( see Ref. [5]) ; 

5. In general, everything that we know about SC is included in Eq. (16) ; 

6. As has been mentioned, Eq. (16) was reproduced by the effective Lagrangian 
approach (see A. Kovner talk at this conference [10] . 

9. The theory status of the equation 

Here, we list all small ( large ) parameters that have been used to obtain Eq. (16): 

as Hl/x) HQ^Ql) ^ 1 ; (17) 
as Hl/x) < 1 ; as HQVQI) < 1 ; (18) 
as < 1 . (19) 

• These equations (Eq. (17) - Eq. (19) ) give the conditions when we can absorb 
all {as ln(l/a;) \n(Q'^/Ql) - terms in the gluon structure function both in 
Eq. (16) and in Eq. (12). 



1 

K < 

as 



(20) 



• This is the main hmitation of our resummation which stems from the contri- 
bution of so called enhanced diagrams ( see Ref. [1] for detail). 

A^c > 1 ■ (21) 

• where Nc is the number of colours. We follow the Veneziano principle [12] of 
large Nc resummation: we select the diagrams in leading order with respect to 
( as Nc )" in each topologic configuration (see Ref. [1] and references therein 
for a detail discussion of this point ). 



SOLUTION 

10. Asymptotic solution 

A first observation is the fact that Eq. (16) has a solution which depends only 
on y. Indeed, one can check that k — Kasympiv) is the solution of the following 
equation 

— -'^ \ '^asymp ) ■ \^^) 

The solution to the above equation is 

.sympiy) 



' Kasymp(y=yo) F {n') 

one can see that Eq. (23) leads to 



= y - yo ■ (23) 



i^asympiy) ^ ^ ln(Q;sy) at asy > 1; (24) 
or xG{x,Q^) oc R^\n{l/x) \n{asln{l/x)) at x — > 

We claim this solution is the asymptotic solution of Eq. (16). To prove this we 
have to consider the stability of the asymptotic solution. It means, that we solve 
our general equation looking for the solution of the form 

1^ = Hasympiy) + Afi;(|/, ^ - ^o) , (25) 

where Ak, is small {An -C Hasymp) but an arbitrary function at { = ^o- We have 
to prove that Eq. (25) will not lead to big An {Ak ^ i^asymp) at large ^. The 
following linear equation can be written for An 

d'AK{y,0 dAK{y,0 _ dF{K) 



In Ref. [1] it was proven, that the solution of Eq. (25) is much smaller than k, since 
^ ^ at large y. 

Therefore, the asymptotic solution has a chance to be the solution of Eq. (16) 
in the region of very small x. To prove this, we need to solve the equation in the 
wide kinematic region starting with the initial condition. We managed to do this 
only in the semiclassical approach. 

11. The semiclassical approach 

In the semiclassical approach the solution of Eq. (16) is supposed to be of the 
form 

K = , (27) 
where S* is a function with partial derivatives: || = a; and f| = 7 which are 



smooth functions of y and ^. 

For an equation in the form F{^,y, S,j,u!) = 0, we can introduce the set of 
characteristic lines {^{y), S{y),uj{y),^{y)), which satisfy a set of well defined equa- 
tions (see, for example, Refs. [13], [14] for the method). In Ref. [1] we developed a 
detailed calculation of the solution. The set of equations is: 

^ <!>{S) . dS ^ 27 + 1 . . 

dy (7 + 1)'' dy {^ + l)^^^'dy 7 + I' ^ ^ 



where $5 = -g^ . The initial conditions for this set of equations 



dlnKin{yo,0 



So = luKiniy 0,^0) ; 7o = 



(29) 



The main properties of these equations have been considered in Ref. [1] analyti- 
cally. Here, however, we restrict ourselves mostly to the numerical solution of these 
equations. 

Fig. 4 shows the behaviour of trajectories and 7 value for the solution. One can 
see that the trajectories can be divided in two groups: (i) first with 7 < — i, and all 
of them approach the trajectories of the DGLAP evolution equation at large Q'; 
and (ii) second one with 7 > — | and the solution of Eq. (16) on these trajectories 
tends to the asymptotic solution at y ^ 1. 

We consider this figure as a nice illustration of our statement that the high energy 
(low x) limit of the general solution to Eq. (16) is the asymptotic solution. 

12. Comparison with other approaches 

In Fig. 5 we plot our calculation for the DGLAP evolution equation, for the 
Glauber - Mueller approach, for the GLR equation and for the new evolution equa- 
tion. We can conclude, that: 



1. SC even in the Glauber - Mueller approach are essential for the gluon structure 
function ; 




FIGURE 4. The trajectories (a) and the 7 values (b) for the solution of the generalized evolution 
equation for nucleon with E? = 5GeV~^. 




FIGURE 5. The current situation in the battle for hdQCD. 



2. The Glauber - Mueller approach considerably underestimates the value of SC ; 

3. The GLR equation leads to stronger SC than the solution to Eq. (16) ; 

4. The new evolution equation does not reproduce the saturation of the gluon 
density in the region of small x which the GLR equation leads to . 



SUMMARY 



13. Where are SC? 

We firmly believe that the new evolution equation gives the correct way of eval- 
uation of the value of SC. However, the difficult question arises: why the SC have 
not been seen at HERA? Our answer is: 

1. The value of SC is rather large but only in the gluon structure function while 
their contribution to F2 is rather small [1], and cannot be seen on the background 
of the experimental errors ; 

2. The theoretical determination of the gluon structure function is not very 
presize and we evaluate the errors as 50% [15] and the SC in xG is hidden in such 
large errors ; 

3. The statement that SC have not been seen is also not quite correct. Our 

estimates [16] show that the contribution of SC is rather large in the F2 slope 
and incorporating the SC one can describe the recent experimental data [17] on 
(see Ref. [16] for detail ). 
14- Last words 

I hope that I gave you a brief but to some extent full review of our approach. 
Obviously, you can obtain more information from two large papers of Ref. [1]. 
Deeply in my heart, I firmly believe that there will be a bright future for the effective 
Lagrangian approach in which one combines the physics of hdQCD, with the formal 
methods of the quantum field theory, and gives a way to incorporate the lattice 
calculation for the non-perturbative observables. However, the parameter or better 
to say a newscalc Ql{x) which is the solution of the equation k,{x,Qq{x)) = 1 
only can be found in the pQCD motivated approach, since the effective Lagrangian 
was only derived and justified assuming this new scale. The fact that the correct 
equation was first found in this approach is also not accidental . since the pQCD 
motivated approach has a great advantage : clear understanding of physics of 
hdQCD. 
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